Ground state configurations of vortices in superconducting film with magnetic dot 
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We consider a thin superconducting film with a magnetic dot with permanent magnetization (nor- 
mal to the film) placed on it by a method based on London-Maxwell equations. For sufficiently high 
dot magnetization a single vortex appears in the ground state. Further increase of magnetization 
is accompanied with the appearance of antivortices and more vortices in the film. We study ana- 
lytically conditions for the appearance of a vortex-antivortex pair for a range of parameters. The 
phase diagram with diversity of vortex-antivortex states is calculated numerically. When appear in 
the ground state, antivortices are at distances comparable to the dot radius. For not too large dot 
radii the total vorticity in the ground state is predominantly zero or one. Magnetic field due to the 
dot and vortices everywhere in space is calculated analytically. 

PACS numbers: 74.81.-g,75.70.-i,74.25.Dw,74.78.-w 



o 

o 
i 

Oh! 

•4-* 



C 

o 

(N 
> 

in 

On 
00 

o 



X 



I. INTRODUCTION 

Type-II superconductors in an external homogeneous 
magnetic field are much studied and well understood*^ 
They accommodate regularly distributed flux tubes- 
vortices for fields between the lower and the upper critical 
fields Similar situation occurs in type-II superconduct- 
ing films in perpendicular fields.—^ A difference is that 
the first critical field for films vanishes for macroscopi- 
cally large samples £JL 

An interesting class of systems which have at- 
tracted a great attention only recently are ferromagnet- 
superconductor hybrid systems. There one examines the 
influence of a material with heterogeneous magnetization 
on a superconductor X Direct contact between the mag- 
netic material and the superconductor is avoided usually 
by a thin insulating layer which suppresses the proximity 
effect. These hybrid systems can be fabricated fully con- 
trolling their parameters^ Inhomogeneous magnetization 
of ferromagnets generates magnetic field that penetrates 
the superconductor. As a response to that field, super- 
currents (and vortices under certain conditions) are in- 
duced in the film. The magnetic field of supercurrents 
interacts with the magnetic subsystem. Therefore by 
tuning parameters of the magnetic subsystem we exam- 
ine different phenomena of the composite system. Hy- 
brid systems offer a number of realizations of new in- 
teresting phenomena which include pinning of magne- 
tization induced vortices, commensurability effects (be- 
tween external magnetic field and periodic structure of 
magnetic material) on the superconductor resistivity and 
others£a£l& 

In this paper we are focused on the simplest 
ferromagnet-superconductor hybrid system which con- 
sists of a single magnetic dot grown on top of a type-II 
superconducting film. The magnetic dot is assumed to 
have permanent magnetization normal to the film sur- 
face. Despite its simplicity this system deserves atten- 
tion since a diversity of different ground states can be 
realized in the parameter space. Dots with sufficiently 
small magnetization produce only supercurrents in the 



film, while with the increase of magnetization differ- 
ent configurations of vortex states appear in the ground 
state(GS)Ji±2i±&iii±^ This is in striking contrast with 
respect to films in homogeneous fields^^ However, we 
expect that for large dot radii any magnetization induces 
vortices in the film, like the case with homogeneous mag- 
netic field. 

The main goal in this paper is to define regions in 
the parameter space with different ground state config- 
urations of vortices. A possible question that is under 
debate in literature is whether antivortices may also be 
induced in the film. One finds different statements about 
the presence of antivortices: while some authors claim 
its existence ; 12 ' 1 ^ 17 the others do not find itj 13 ' 14 We 
think that rough estimates in Ref^, calculations with 
magnetic dipolc in RefJ^ or study inside the nonlinear 
Ginzburg-Landau theory with restriction to zero total 
vorticity states done in RefJ£ are insufficient, and offer 
an independent and detailed study of this problem. Un- 
der which conditions and where vortices and antivortices 
appear will be answered in our paper. 

Magnetization of the dot (which is normal to the film 
surface) produces magnetic field in the film, which is un- 
der the dot parallel to the magnetization, but antiparallel 
for larger distances than the dot radius. The magnetic 
field under the dot favors the appearance of vortices un- 
der the dot, the appearance of antivortices outside the 
dot. Whether some vortex-antivortex configuration have 
the smallest energy or not depends on details of the in- 
teraction energy. 

Vortices have spatial structure and the superconduct- 
ing order parameter vanishes roughly at distances smaller 
than £ around the center of the vortex. Here £ denotes 
the coherence length of a superconductor X The nonlin- 
ear Ginzburg-Landau approach takes this into account, 
but apart from numerical treatment one can hardly get 
analytic results^ Our approach is based on London- 
Maxwell equationsi 13 ' 18 Despite its simplicity (it treats 
vortices as point objects) it an useful approach since one 
may get analytic expressions for relevant quantities which 
are asymptotically exact when all lengths in the problem 
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are larger than the spatial vortex extension. 

It is interesting to mention that the vortex nucleation 
in superconducting microtriangles and squares occurs in 
such a way that the symmetry of the superconductor is 
preserved. 19 ! 20 For example in the case of a triangle, the 
state with total vorticity two is realized with three vor- 
tices and a single antivortex in the center of the triangle 
and preserves C3 symmetry. Our system with infinite 
superconducting film and with magnetic dot on top of it 
has rotational symmetry around the dot center. This 
symmetry is reduced in the presence of vortices and an- 
tivortices contrary to the abovemcntioned example. For 
the most possible symmetric states with total vorticity 
zero the rotational symmetry is either reduced to dis- 
crete Cn symmetry (for the case of N > 2 single vortices 
and antivortices) or does not exist at all (for a vortex- 
antivortex pair). When two and more (we have checked 
up to four) vortices and antivortices are present in the 
zero vorticity state, vortices are distributed in a sym- 
metric fashion around the dot center while antivortices 
are outside the dot in the same manner, like homotheti- 
cally transformed vortices. In the case of a single vortex- 
antivortex pair, the dot's center, vortex and antivortex 
are collinear. Other states with nonzero vorticity are less 
symmetric. When the number of antivortices is five and 
more, they may form shells around a central vortex in 
zero vorticity states*^ 

The rest of the paper is organized as follows: in section 
II we introduce a theoretical model we use to describe 
the system. We calculate analytically the interaction en- 
ergy between a cylindrical magnetic dot and vortices in 
the film. In section III we determine analytically con- 
ditions for the ground state configurations that consist 
of a single vortex and a vortex-antivortex pair, as well 
as separatrices between those ground states. Positions 
of vortices are also found. Magnetic field in whole space 
and supercurrents in the film are calculated in section 
IV. Numerical results for the phase diagram with diver- 
sity of vortex-antivortex configurations as well as their 
positions are presented in Section V. Section VI contains 
numerical estimates of the parameters and conclusions. 



Some technical details are relegated for the appendix. 



II. MODEL 

We consider a circular magnetic dot of radius R and 
thickness at placed above the infinite superconducting 
film at distance d with its basis parallel to the film sur- 
face. The dot magnetization M is assumed to be con- 
stant and normal to the film surface, see Fig. [T] Appart 
from the supercurrents, the dot may also induce vortices 
and antivortices in the film. We will study these vortex 
configurations. Since our problem has many parameters 
there are many regions in the parameter space which 
may have different ground state configurations(GSC). 
The GSC here denotes a configuration of vortices in the 
film with lowest energy. Among all possible GSCs a triv- 
ial one has no vortices. This is expected for sufficiently 
small magnetization. With increasing M the appearance 
of vortices is energetically favorable. 

We assume quite generally that our system consists of 
N vortices with vorticities rii at positions p i for a given 
magnetization M. The energy of the system is 

N N N 

E ni ,...n N = y~]njU v + y^^n i njU vv (pjj) + ^ njU mv (pi), 

i—l i<j i—1 

(1) 

where pij = \p i — Pj\ and p i = \p i \. U v is the single vor- 
tex energy, U vv is the vortex-vortex interaction, and U mv 
is the vortex-magnet interaction. For a given magneti- 
zation of the dot, the system will be in the state where 
E ni .... nN is minimal. Expression ([1} assumes that the 
system in the trivial ground state has zero energy Eq = 0. 

U mv and U vv can be calculated by using the approach 
developed i n 13 i 18 based on London-Maxwell equations. 
Here we only quote final results. For the interaction en- 
ergy between the dot and a vortex with vorticity n placed 
at distance p from the dot's center we get 



U mv (p,n) = -nMR(j>o 
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where Jo and J\ are the Bessel functions of the first kind. A is the effective penetration depth and is equal X 2 L /d s , 
where Xl is the London penetration depth and d s the film thickness. The film thickness is assumed to satisfy d s <C Xl 
and accordingly the magnetic field and currents are uniform through the thickness of the film. Expression @ has the 
following behavior in some regions (we take d = for simplicity) : 
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FIG. 1: Magnetic dot with perpendicular permanent magne- 
tization M upon infinite superconducting film. 
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Some technical details for obtaining ^ are presented in 
the appendix. 

The interaction between vortices of vorticities n\ and 
n.2 separated by a distance p is given by 



U vv (p) 



nin 2 <f>o 



16?rA 

which has asymptotic forms 

2U V , p<£<A 

fv,(//i ■nn 2 { Sx( A -T ')■ - ■ I''' A (5) 

A < p 
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In previous expressions 0o is the magnetic flux quan- 
tum. H is the Struve function of order zero and Y is 
the Bessel function of the second kind of order zero.— U v 
is the single vortex energy 2 ^ 



U v = 



(6) 



where 7 « 0.577 is the Euler constant. Two vortices 
which centers are at distances smaller than £ are consid- 
ered in our model as a double vortex. Its total energy is 
4U V : from two single vortices comes 2U V as well as from 
their interaction. In the same way a vortex-antivortex 
pair at distances smaller than £ does not exist, i.e. it is 
annihilated. Here we mention that Eqs. (J4j) and ((6|) are 
valid for films with lateral dimensions much larger than 
the effective penetration depth. In the opposite case the 
lateral system size enters expressions ((4]) and ([6]) instead 
of A. This change does not complicate further consider- 
ations and we do not consider such situation. 

In the following we will use energy expression ([T|) which 
will be minimized in the parameter space and which de- 
termines the structure of vortex configurations in the GS. 



III. GROUND STATES WITH VORTICES 

To get further insight into possible GSCs of the system 
with respect to different parameters, we consider asymp- 
totic forms of Eqs. {2]) and (j4|). To simplify expressions we 
will just consider the case of thin magnetic dots (at < R) 



FIG. 2: (Color online) Phase diagram of a magnetic dot for 
L — 6 and a t /R = 0.01: vortex-antivortex state appears 
for large enough dot magnetization and large enough X/R. 
Solid lines are obtained numerically, while dashed are ana- 
lytic formulae 0, @ and JTUJ with A = 1.7. Inset: up- 
per(lower) curves show antivortex(vortex) positions. Solid 
(Dashed) lines correspond to the upper (lower) phase bound- 
ary of the vortex-antivortex pair. Dotted line are plotted 
using analytic formulae (|12[) and (I13|l . 



placed at the film surface (d = 0). The other cases may 
be straightforwardly done having the asymptotic forms 
of U mv . 

Necessary condition for the appearance of an extra vor- 
tex with respect to the trivial GSC is E\ < Eq. Using 
© in the lowest order in A/ R we get 



. , 2A 



(7) 



Here we have introduced M. = Mat(f>o/U v . This result 
agrees with one from Ref^. Here we have taken into ac- 
count that the strongest attraction energy between single 
vortex and magnet occurs at p = 0. The phase bound- 
ary is a linear function of X/R for R < X. For large dot 
radii R, result ([7]) shows that any nonzero magnetization 
induces a vortex in the film. This resembles the result 
for vanishing of the first critical field for thin films in an 
external homogeneous magnetic field^^ 

Apart from the trivial GSC and the single vortex state, 
there are other states for higher M.. We will now deter- 
mine a portion of the parameter space where a vortex- 
antivortex pair appears in the GS. We found that such 
possibility is only briefly mentioned in literature ; 12 ' 17 or 
disproved^ We will try to clarify this issue giving ex- 
plicit expressions for the phase boundary and for the 
vortex and antivortex positions. 

The energy 2?i _i of a vortex-antivortex pair separated 
by P (p <C 2A) with vortex under the dot center reads 



M 



AXp 



E h -l{p) Umv(P,l) 

u v u v 

2 -R 2 \ / 4A ^ 
L +M 8V) [ ln J-l 



(8) 
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where L is the logarithmic factor in the vortex self en- 
ergy L = ln(4A/£) — 7. We may notice from Eq. 
that the main contribution from the magnetic dots on 
the antivortex is the energy cost which scales with the 
distance from the dot as 1/p, while the energy gain due 
to the vortex-antivortex attraction scales logarithmically 
which may lead to a stable potential minimum for some 
parameters. For that energy minimum we get that it 
occurs at p* = 2A/(1 + j- ^li 1 ) wruc h should be sup- 
plemented with the self consistency condition p* <C A 
where Eq. ([5]) is valid and also p* > a t , R. Implicit 
equation Ei t _i(p*) = defines the phase boundary for 
the creation of a vortex-antivortex pair. In addition if 
< Ei and £4.-1 < are satisfied the vortex- 
antivortex pair forms the GS. Condition E\ _i < E± gives 



M < 



X 2 



32 



R 2 L [exp (l + 7 +f)-2] 



(9) 



which is the upper dashed line for vortex-antivortex 
phase boundary in Fig. [2] Another phase boundary we 
get from the condition that the antivortex is outside the 
dot p* > AR which gives 



, , 8 A 
M> -LR A 



(10) 



with some number A of order one. 

From the last two inequalities we get a condition on 
X/R when the vortex-antivortex forms the GS (in the 
case a t < R, d = 0): 



\_ A 

R ~ 4 



<'xs>(l + 7 +-] 2 



(11) 



In Eq. ([5]) we have assumed that the vortex posi- 
tion is under the dot pi = for simplicity. Having 
in mind relatively flat magnet-vortex interaction ([3]) for 
p < R the system can gain even more energy allowing 
pi > toward antivortex, since the energy loss in the 
vortex-magnet interaction may be overcompensated by 
the vortex-antivortex interaction . A similar calculation 
to the previous one gives for the vortex displacement 



pi _ 64A 2 
~R ~ M 2 R 2 L 2 

while for the antivortex 
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and the center of the dot, vortex and antivortex are 
collinear. Physically we see that stronger the magnetiza- 
tion of the dot. the vortex is closer to dot's center, which 



is expected. What we also see is that the symmetry of a 
single vortex GSC is violated for the range of parameters 
where the antivortex appears. 

IV. MAGNETIC FIELD 

Magnetic dot on top of a superconducting film induces 
circular supercurrents in the film. These currents gener- 
ate magnetic field in and outside the film. Total magnetic 
field in space is a sum of three terms: due to supercur- 
rents, due to the dot and due to vortices. In this section 
we calculate the magnetic field using the approach de- 
veloped i n 13 i 18 . The presence of vortices in the film may 
be observed by measuring the magnetic field and its be- 
havior near the film surface, since vortices change the 
magnetic field dependence of radial separation from the 
dot's center. 

A vortex of vorticity n produces normal to the film 
surface (axial) and parallel the film surface (radial) mag- 
netic field which respectively react^ 3 . 



k exp(— k\z\) 



U v f \ n< Po f ,, fees 
B Z (P,Z) = —J Q dk — 

BHp, Z )= n ^r^ 

11 2?r z J 1 



2Afc 
k exp(— k\z\) 



2Xk 



Jo(kp), (14) 
Ji(kp). (15) 



The previous expressions at the film surface z — have 
the following asymptotic forms: 



BUp,z 




(16) 
(17) 



The parallel magnetic field outside the film (remember 
that in our calculations the film is just in z = plane) 
changes the sign going from one side of the film to an- 
other. This is due to the fact that the vortex induces 
currents in the film which circulate around it, and jump 
in BV crossing the film surface is the condition for the 
jump at the boundaries in electrodynamics due to surface 
currents^ Normal magnetic field is continuous across 
the film. The surface current density K v produced by 
the vortex is given by 

K„ = ^-e z x (B|(p, 0+) - Bf (p, 0-)), (18) 

where c is the velocity of light. 

The magnetic field due to the dot is given by two in- 
tegrals: 
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B*(p,z) = -1%MR j dkJ (kp)J 1 (kR)^^j 



exp(-k(\z\ + d)) 



2Xk 



[1 — exp(— ha)} 



Bf(p,z) = -2irMR^ d/cJi(fcp)Ji(/ci?){sign(2)^Y 



exp(-k(\z\ + d)) 



2\k 



1 — exp(— ka)] 



+ cxp(— k\d — z\) — exp(— k\d + a — z||. 



(19) 



+ sign(d — z) [1 — exp(— k\d — z\)] — sign(d + a — z) [1 — exp(— k\d + a — z\)] i, 



(20) 



Eq. CL9])(Eq. d20j) ) we write as a sum Bt^ = B 



of fields due to the dot B 



din 



dm 

KID" 



z(ll) 



id dm 

B <\\) 



and due to the supercurrents 

^(11)' * s f° rmau y defined by setting A — » oo in 

Eq. (Un])(Eq. (|2H)l ) which means the system with the dot 
and without the film. We evaluate (fT9|) and (f20|) for d = 
and at the film surface. Purely magnetic terms are given 
by 



B d z m {p 7 0) = ~TrMR 2 a t ^r, 



B dm (p,0) = -nMR 2 a t ^, 



(21) 
(22) 



to the leading order for p » a tl R. The previous result 
can be understood as the magnetic dipolar field. 24 The 
part due to supercurrents is 



Sf(p,0) 



<(p,z^0) 



-irMR a t 



Again parallel field B^ jumps across the film surface due 
to surface currents. The surface current density K m in 
the film due to the presence of the dot is given by 



a t ,R< /?< A 
a t , J R< A < p 

(23) 

jr, a t ,i?< i o<A 
S a t ,i?<A<p 
(24) 



K m =-^e z x (Bf(p,0+) 
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Bf(p,0")). 



(25) 



As first proposed in Refj 1 ^ the presence of a single vor- 
tex in the film can be proved by observing the change of 
sign of the total field B z near the film. We can now cal- 
culate that it happens when B z (p sz ,0) + B z (p sz ,0) = 
or at p sz — 2iry/ M R 2 a t X/4>o provided a single vortex ap- 
pears in the GS. Using the condition for the single vortex 
(0 we get p sz « ^/LRX/2. We can also formulate a sim- 
ilar condition for the presence of a single vortex but for 
the parallel field, which also changes sign for at distance 
p s || from the dot defined by B$(p s pO) + BV,(p s pO) = 0. 
The solution of the resulting cubic equation is p s \\ « 

(67T 2 Mi? 2 a 2 A/0 o ) 1/3 above the film (z -> 0+), which us- 
ing becomes p s \\ w (3LRa t A/4) 1 / 3 . 



B,(p,0)~-l/p 3 

dipolar field 
—\ h 



l/p 
vortex field 



1/P 3 
vortex field 



at, R 
Bn(p,0) 



+ 



-i/V 

dipolar field 



Up 

vortex field 
H H 



Up" 

vortex field 



at, R 



P4 



X 



FIG. 3: Behavior of the magnetic field near the upper film 
surface. Magnetization of the dot is assumed to be directed 
as in Fig. [T]such that a single vortex appears under the dot. 
The magnetic field changes the sign at at some distance from 
the dot due to the presence of the vortex. 



Qualitatively different behavior of the normal and par- 
allel magnetic field is summarized in Fig. [3] when a single 
vortex is present in the GS. Sufficiently close to the dot 
dominates the dipolar field from the dot, while at larger 
distances the vortex part of the field is a leading term. 

So far we have only considered a single vortex un- 
der the dot. Since there may be many other vortex- 
antivortex states in the film, the magnetic field (vor- 
tex part only) of such configurations behaves differently 
than in (fTT))) and ([IT)) and will get angular dependence. 
We do not analyze this case. Again, the total magnetic 
field close enough to the dot will be dominated by the 
dipolar field from the dot, while at larger distances the 
anisotropic vortex part of the magnetic field prevails. 
One should be able to use the measurement of magnetic 
field near the film for detection of many-vortex states 
bounded by the dot. 



V. NUMERICAL STUDY OF GROUND STATES 
WITH LOW NUMBER OF VORTICES 

In section III we have shown analytically a possibility 
of having a vortex-antivortex pair in the GS. That was 
the simplest GSC with antivortex. Certainly there are 
other more complicated states for a range of parameters. 
Analytic study of these states is in principle straightfor- 
ward using already introduced expressions, but tedious. 
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FIG. 4: (Color online) Phase diagram of a magnetic dot for 
L — 6 and at/R = 0.01 with different vortex-antivortex con- 
figurations. 
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FIG. 5: (Color online) Phase diagram of a magnetic dot for 
L = 6 and a t /R = 1 with different vortex-antivortex config- 
urations. 



In this section we study numerically GSCs. For a given 
configurations uniquely determined by m,7i2,.. ■ we use 
expression ([I]) to obtain the energy minimum and po- 
sitions of vortices and antivortices. The GSC for given 
X/R, at/R, L and M. has the energy minimum over of all 
possible vortex-antivortex configurations. In our calcu- 
lations we took into account just states with low numbers 
of vortices, since they cover significant part of the param- 
eter space as well as the other states are computationally 
demanding. Dot-magnet distance is set to d — 0. We 
have examined all states which have up to four vortices 
with antivortices, while for the states with five vortices 
we took into account vortices without antivortices. This 
is certainly not correct since antivortices appears in states 
with five vortices as well, but such states are located in 
a particular region of the parameter space and do not 
affect states with up to three vortices (see further in the 
text). 

States with six and more vortices are not taken into 
account. This is not crucial for our study since the vor- 
tex state with vorticity n will appear in the phase dia- 
gram for M. > 2n\/R which is either for large magneti- 
zation or large magnetic dots (with respect to A) . On the 
other hand line with n single vortices (and antivortices) 
is moved a little bit toward smaller M. for a given A/i?. 

In Fig. [4] we show the GSC of the magnetic dot with 
a t /R = 0.01 and L — 6. For low M. there are no vor- 
tices in the GS for any X/R. With increasing M. state 
Ei appears. Further increase of M. is ultimately followed 
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FIG. 6: (Color online) Vortex distances (in units of dot radius 
R) from the dot center for L = 6 and at/R — 0.01 which 
correspond to a particular ground state configuration given 
in FIG. H 



with the antivortex appearance in the GS. To be concrete 
let us consider X/R = 50. After E\, state with antivor- 
tex Ei t -i appears which is the GSC for some region of 
M. Then states £a,i,-i,_i, ^i, 1,1,-1,-1,-1, ^1,1,1,-1, -ij 
-Ei, 1,1, 1,-1, -1,-1,-1 appear for larger M. Further states 
might have included states with five vortices with an- 
tivortices if we had taken them into account. We see that 
states that evolve from each other have vortex and/or 
antivortex more/less with respect to its neighbors. We 
believe that this should mean that states with three vor- 
tices will not be affected by the states with five vortices. 

Also the net vorticity of the GSCs is in general small 
(0,±1), just near the origin it is higher. However for 
small X/R the dot radius is pretty large for realistic films 
and further analysis with finer resolution of X/R is nec- 
essary. We are here mainly concentrated in the region 
where R < X since it is the most interesting experimen- 
tally, see the last section for some numerical values of 
parameters. 

The distance between the center of the dot and vortices 
and antivortices are shown in Fig. [6] and Fig. [7] respec- 
tively. We see that they abruptly change at boundaries 
between different GSCs. We also see that antivortices 
when are present are at distances of a few R. Inside the 
same GSC for fixed X/R with increasing M. the antivor- 
tices spread, while vortices shrink. This is plausible: very 
large magnetization of the dot would expel antivortices to 
very large distances since the interaction energy cannot 
be then compensated by the attractive vortex-antivortex 
attraction. On the other hand vortices are attracted at 
smaller distances toward the potential minimum of the 
vortex-magnet interaction. These results are also con- 
firmed by analytic formulae (fT2|) and (fT5|) from the study 
of a single pair. 

We also find that in states with zero total vorticity and 
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FIG. 7: (Color online) Antivortex distances (in units of dot 
radius R) from the dot center for L — 6 and at/R = 0.01 
which correspond to a particular ground state configuration 
given in FIG. Q] 
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FIG. 8: (Color online) Relative energy change in the ground 
state when antivortices are included for L — 6 and at/R = 
0.01. 



two, three and four vortices, positions of antivortices are 
obtained as homothetically transformed positions of vor- 
tices, and they form line, equilateral triangle and square, 
respectively. 

The GSCs without antivortices are shown in Fig. [9] 
Relative energy difference due to the presence of antivor- 
tices is shown in Fig. [8] The relative energy gain is quite 
significant which means that states without and with an- 
tivortices have quite different energies. 

We see from Fig. [8] that for smaller magnetization the 
vortex states with single vortices appear, while state with 
double vortex Ei is present for higher M. Simple ana- 
lytical check gives that condition E\ t \ < E% translates 
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FIG. 9: (Color online) Phase diagram of a magnetic dot for 
L — 6 and at/R — 0.01 without antivortices. 



into 



M < - exp (2L + 1 - 27) 



R 

A' 



(26) 



which for L = 6 becomes M < 23200i?/A. That also 
gives one explanation why the only ground state with 
giant vortices in Fig. [4] is E2,—i,—i' it appears around 
X/R = 100, M = 250 in agreement with rough estimate 
(|26|) . The other states with giant vortices may appear 
only for higher M.X/R ratios. 

For comparison with thicker magnetic dots we have 
calculated the GS for at/R = 1 and they are shown in 
Fig. O The single vortex appears now for larger A4, 
which is obvious from ([7]). State £1,-1 now occurs for 
larger X/R. However a global picture with diversity of 
vortex-antivortex states again holds. However GSC here 
are shifted toward higher X/R and M. values with re- 
spect to the corresponding states for thinner dots. We 
may interpret this as a fact that thicker magnets are less 
efficient in exciting vortices which is due to the larger 
extent of the dot with respect to the film surface. Very 
thick dots [at 3> R) interact with vortices weakly, see 
two cases in formula ([3]). 

Now we comment about the applicability of the Lon- 
don approximation for our case. The dot radius over 
the coherence length £ as a function X/R is given as 
R/£ = Rcxp(L + 7)/(4A). We may conclude that for 
X/R < 100 and L = 6 the dot radius is always larger 
than the coherence length. On the other hand a neces- 
sary condition for the London approach to be valid is that 
all lengths are (much) larger than the coherence length. 
We see that this is satisfied better as X/R approaches 
smaller values. We expect that our results with different 
GSCs are valid even qualitatively for X/R < 50, while the 
general picture with antivortices and low vorticity GSC 
holds even for smaller dot radii. However it is expected 
that for smaller R/^ ratios, giant vortices are favorable. 
We have already mentioned such tendency. 

Let us mention that a certainly better numerical way 
would be the nonlinear Ginzburg-Landau theory. For the 
case of mesoscopic superconducting discs comparison of 
these two methods shows that both methods give similar 
results for R = 6£ and up to five vortices^ 
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VI. DISCUSSIONS AND CONCLUSIONS 

Let us consider numerical values of parameters for re- 
alistic systems. For thin Nb films close to the critical 
temperature T/T c = 0.98 values for the London penetra- 
tion depth and the coherence length are A^ = 560 nm 
and £ = 58 nmS Then L = 6 corresponds to a film 
of thickness d s = 30 nm with the effective penetration 
depth A = f 0.5 jum. Condition for state (TTTT) can be 

rewritten as R < Al y£7 (9^s) or in our case R < 260 nm. 
The condition for the applicability of the London theory 
£ < R is still satisfied, so the results should be valid. 
Larger magnetic dots may have other vortex-antivortex 
configurations, and the London theory is expected to be 
applicable. 

In this paper we have considered infinite films. We 
expect this not to be a severe limitation as soon as the 
dot radius is much smaller than the system size since the 
antivortices when appear are at distances of the order of 
R and the boundaries of the film should not affect it. For 
films with lateral dimensions smaller than A the single 
vortex energy, Eq. ([6]) will have under the logarithm the 
lateral dimension instead of 4A, and for small dot radii 
essentially the same story should be repeated, just with 
new U v . 

To conclude, in this paper we have considered a 
thin superconducting film with a cylindrical magnetic 
dot with permanent magnetization upon it. Inside 
the Maxwell-London approach we have calculated the 
vortex-magnet interaction as well its asymptotic limits. 
Using these results we have shown analytically that a 
vortex-antivortex pair appears in the ground state of 
the system for some range of dot's radii and its mag- 
netization. Necessary magnetization for that is compa- 
rable to the magnetization for the appearance of a single 
vortex. Magnetic field everywhere in space is also cal- 
culated. Near the film surface, the magnetic field has 
different scaling forms with the distance from the dot 
center. This fact may be used for the experimental de- 
tection of vortices. In addition to that we have calculated 
numerically the phase diagram with up to four vortices 
with antivortices. 
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VII. APPENDIX 

In this appendix we calculate asymptotically the in- 
tegrals that appears in the interaction energy We 
consider an integral of the form 

exp(— cx) 



I (a, b, c) 



dxMax)Mx) x{1 + 2bxy (27) 



In the limit a 3> c and a 1 this integral will be cut 
by the oscillations of the Bessel functions and the main 
contribution comes from the region x < 1/a <C f. Then 
we can expand Ji(x) exp(— cx) « x(l — cx)/2. Using the 
tabulated integrals^ 



1 J (ax) it vim 
dx 1 + x = 2 t H °( a ) ~ 



,oo x 

I dxJo(ax) = — , 
Jo " 

we easily get 



(28) 
(29) 



I(a,b,c) 



2b 



8b 2 2 
Using the expansion^! 

I [Y (x) - H (x)] -- 



Ml) "Ml) 



55- (30) 



7 + log f- x, x<l , . 
-x^+x- 3 , x>l 



where 7 rs 0.577 is the Euler constant, one can further 
simplify the asymptotic expressions for /(a, 6, c). 

In the limit a <C c and a <§; 1 integral (f2"T| is cut by 
the exponential function at x ~ 1/c, which means the 
argument of Jo function ax ~ a/c -C I and we expand 
J(0,ax) w 1 — a 2 x 2 /A. We get 



I(a,b,c) = J dxJi(x) exp(— cx) 
b x l + 2bx 



(32) 



Then using the tabulated integral 



,27 



dx^ — '- = 1 + - + - [Yx(a) - ffi(a)] 



1 + x 



dxJi (ax 



CX P (-7X) r— 

) = -7+ Vl + 77 



poo 

/ dxJi(ax) exp(— jx) = 1 -= 

Jo ' \/l 



and the expansion for 



l[H 1 ix)-Y 1 (x)] = ~ + Z{ 1 -2'y + ] g± 



we get for a -c c <C 1 <C b 
1 



I (a, b, c) 



1 - 27 + c(a 2 - 4) 



V 1 + c 2 + log + c- 



(33) 
(34) 
(35) 

(36) 
(37) 



I the above expressions Jo and J\ are the Bessel functions 
of the first kind, Yq and Y\ are the Bessel functions of the 
second kind, while Hq and Hi are the Struve functions 
of order zero and one, respectively^ 
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